Transport signatures of spatially modulated electronic nematic phases 



Hae- Young Kee/' ^'0 Christoph M. Puetter,^ and David Stroud"^ 

^Department of Physics, University of Toronto, Toronto, Ontario MBS 1A7 Canada 
^Canadian Institute for Advanced Research, Toronto, Ontario Canada 
^ Craduate School of Pure and Applied Sciences, University of Tsukuba, Tsukuba, Ibaraki 305-8571, Japan 
^Department of Physics, Ohio State University, Columbus, OH U.S. A 

Electronic nematic phases are broadly characterized by spontaneously broken rotational symme- 
try. Although they have been widely recognized in the context of high temperature cuprates, bilayer 
ruthenates, and iron-based superconductors, the focus so far has been exclusively on the uniform 
nematic phase. Recently, however, it was proposed that on a square lattice a nematic instability in 
the d-wave charge channel could lead to a spatially modulated nematic state, where the modulation 
vector q is determined by the relative location of the Fermi level to the van Hove singularity. [T] 
Interestingly, this finite-q nematic phase has also been identified as an additional leading instabil- 
ity that is as strong as the superconducting instability near the onset of spin density wave order. 
[2] Here we study the electrical conductivity tensor in the modulated nematic phase for a general 
modulation vector. Our results can be used to identify modulated nematic phases in correlated 
materials. 



PACS numbers: 71.30. +h,72.10-d,72.15.-v 

Introduction - Identifying the genuine ground states of 
doped Mott insulators has been one of the key issues in 
the field of strongly correlated electron systems. It was 
proposed that in the ground states of doped Mott in- 
sulators, electrons arrange themselves in certain charge 
and/or spin density patterns due to strong electron inter- 
actions and quantum fluctuations introduced by electron 
or hole doping. [3] These self-organized phases are clas- 
sified by different broken symmetries. Among them is 
the electronic nematic phase, which is characterized by 
a spontaneous Fermi surface (FS) deformation that re- 
duces the rotational symmetry of the underlying lattice. 
On the square lattice the nematic order parameter has 
the form 

An,a = ^(cosA:^ - cosA:^)(4^^Ck,a), (1) 
k 

where c\.^ (ck,cr) creates (annihilates) an electron with 
momentum k and spin cr . For a charge nematic phase one 
has A^^^ = A^^i = A^. Note that a finite represents 
90 degree rotational symmetry breaking, as A^ changes 
sign under such a lattice rotation. 

Combined with the tight binding model for a square 
lattice, a simple nematic mean field Hamiltonian is then 
given by 

Hmf = ^ [ck - /i + An(cos - cos ky)] 4,a^k,a(2) 

k,cr 

where ek = —2t[coskx + cos ky] — At' cos k^ cos ky — 
2t''[cos {2kx) +COS {^ky)] and /i is the chemical potential. 
The parameters t, t', and t" denote nearest, next-nearest, 
and third- nearest neighbor hoppings, respectively. The 
uniform nematic phase can hence be associated with dis- 
tinct nearest-neighbor hopping integrals along x- and y- 
directions. During the last few years, theoretical studies 



have focussed on the microscopic route to such an ef- 
fective Hamiltonian, the experimental consequences, and 
the nature of the isotropic-nematic phase transition. [4- 
15 Fascinating experimental results, which support the 
existence of the electronic nematic phase, [16] were also 
reported in the high temperature cuprates, [TT^^ bi- 
layer ruthenates, [2TH24j and iron-pnictides. ^25-.-.-28j 




FIG. 1: [color online] General nematic modulation (see main 
text for details). 

Although the theoretical focus so far has been lim- 
ited to the uniform nematic phase described by Eq. (1), 
there is no a priori reason why the uniform nematic 
phase should be the leading instability among the possi- 
ble anisotropic instabilities of an isotropic metal. In this 
context, two recent theoretical proposals are interesting 
to note. Metlitski et al. in Ref. [2^ studied instabili- 
ties near the onset of spin density wave order and found 
that modulated nematic bond order is an equally strong 
instability compared to the well-known superconducting 
instability. Holder et al. in Ref. [1 investigated a ne- 
matic instability with a finite wave vector and found that 
the direction of the modulation wave vector depends on 
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the filling relative to the van Hove filling. An extended 
Hubbard model was furthermore studied to analyze its 
renormalization group flow. [29] The spatially modulated 
nematic order suggested in these studies can be described 
by the following bond order parameter operator, 



An(q) 



^(cos 

k,cr 



COS A:2y)cj^^q^2,cr^k-q/2,cr, 



(3) 



where a finite q determines the periodicity and direction 
of the modulation and the form factor (cos kx — cos ky) = 
(ik represents the nematic character. A typical example 
is given in Fig. [l] where the red and blue colors on the 
bonds indicate the nematicity, and the modulation of the 
bonds gives rise to a striped pattern tilted at an angle 6 
relative to the x-axis. 

In this paper, we study the transport behavior of such 
a modulated nematic phase using the effective medium 
theory. [3QH32] Our results can be used to pinpoint the 
existence of the modulated nematic phase in correlated 
materials in general. Before we proceed to the study of 
transport, let us check what factors determine the direc- 
tion of the modulation vector q, as it leads to a qualita- 
tive difference in the resistivities which will be presented 
later. 

Direction of the modulation vector - In Ref. [1 , the 
modulation vector q is determined by the shape of the FS 
close to van Hove filling, assuming that the nematic in- 
teraction has only a weak momentum-dependence around 
this regime of FS. It was found that for fillings above the 
van Hove singularity (vHS) the modulation vector is ori- 
ented diagonal and below parallel to the crystal axes. 
Here we study a tight binding model with a different FS 
shape and find opposite behavior where q is parallel to 
the X- or y-axis above the vHS and lies diagonally below. 

Let us consider an effective Hamiltonian of the form 



k,cr 



(4) 

where ^(q) denotes an effective attractive interaction fa- 
voring a nematic instability and N is the number of lat- 
tice sites. 

Assuming ^(q) does not have a strong q-dependence, 
the nematic susceptibility is determined by the d-wave 
static {uo = 0) bare polarization function (the d-waviness 
originates from d\^) given by 



^ ^2 /(^k+q/2) - /(^k-q/2) ^^-^ 
^ k 



ek+q/2 — ek-q/2 



where /(^k) is the Fermi function and ^k = ^k — Simi- 
larly the s-wave polarization function can be obtained by 
setting dk = 1. 

Setting the tight binding parameters to t = 1, = 
0.375 and t" = 0, Fig. [2] (a) shows the FSs for several 
chemical potentials /i. The values ja = 1.45, 1.5 and 1.52 
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FIG. 2: [color online] (a) Fermi surfaces for different chemical 
potentials /x = 1.45 (blue), 1.5 (black) and 1.52 (red). The 
shaded area indicates the occupied states, (b) The bare d- 
wave (solid lines) polarization function at zero frequency for 
the same /i- values as in (a). The bare s-wave function (dashed 
lines) is also shown for comparison, (c) The bare d-wave po- 
larization function in the two-dimensional q-space for the 
three cases below, at and above van Hove filling with 2/ci?- 
lines (dotted lines) partially superimposed. The 2/cF-lines are 
defined by ^(q+G)/2 = 0, where G belongs to the set of recip- 
rocal lattice vectors. [1 (d) The bare s-wave polarization for 
the same parameters as in (c) for comparison. 



lie below (blue), at (black) and above (red) the vHS, 
respectively. The d-wave (solid lines) and s-wave (dashed 
lines) static polarization functions for each /i along F — 
X — M — r are plotted in (b). In all cases, the d-wave 
polarization is more singular than the s-wave case. At 
van Hove filling, the logarithmic singularity occurs at q = 
leading to the pure nematic phase. However, for ja 
below the vHS (blue solid line), the d-wave polarization 
function is singular at a finite q along the F — M direction 
in q space. On the other hand, above the vHS (red solid 
line), it is most singular for q || T — X. The panels (c) 
and (d) show the d-wave and s-wave bare polarization 
functions in the two-dimensional q-space. 

Qualitatively, the results are the same as those ob- 
tained in Ref. [1 . The main quantitative difference 
though is the relation between electron filling and the 
direction of q. However, this difference is simply due 
to the relative meaning of filled and empty states. For 
electron-filling [the grey area in (a)] the blue FS is be- 
low the van Hove points, but for hole-filling it is above. 
Here as well as in the study fD the preferred orientation 
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of the modulation vector for the blue FS is along the 
diagonal direction. Thus the direction of the modula- 
tion vector depends on the FS shape near the van Hove 
points/saddle points rather than on the position of the 
Fermi level with respect to the vHS. Now let us present 
how the transport in the modulated nematic phase is af- 
fected, in particular, when the leading nematic instability 
occurs for q parallel or diagonal to the crystal axes. 

Conductivity Tensor - The spatially modulated ne- 
matic phase can be viewed as a collection of uniform ne- 
matic states with varying bond strengths and periodicity 
27r/|q|. As mentioned before, an example of the modu- 
lated nematic phase with q = (~|f 7 ff ) shown in Fig. 
[1] where the red/blue bond has a higher/lower hopping 
integral. For a given ^/-position, the red bond changes 
to blue and vice versa along the x-direction. Moving 
from row to row, the modulation along the x-direction is 
shifted by a phase factor. The finite-q nematic phase 
can also be viewed as a collection of parallel stripes of al- 
ternating, approximately uniform nematic domains [33] 
running at an angle 6 relative to the crystal x-axis, since 
the nematic phase possesses Ising symmetry. 

The principal components of the conductivity tensor of 
each nematic domain are denoted by cFxx and (jyy form- 
ing a diagonal 2x2 conductivity tensor. The principal 
axes are in the same direction as the crystal axes, as the 
directions are pinned by the lattice potential. [12 Based 
on the Ising symmetry, there are two types of nematic 
domains. One type has higher Gxx-, while the other has 
lower a XX (and vice versa for cFyy). We denote these two 
types of nematic domains A and B, and assume an area 
fraction of p for type A and 1 — p for type B. Thus the 
conductivity tensor for each domain is given by 



Refs. |30], [3T] and |32]) 



' A/B 



' A/B ^ 



(6) 



where cj^ = ay^ and = ag^. Since the stripes make 
an angle of 6 relative to the x-axis, we introduce the ro- 
tated coordinate system — y' in such a way that the 
x'-direction runs parallel and the ^'-direction runs per- 
pendicular to the stripes (see Fig. [T]). Our goal is to get 
the effective conductivity tensor in the x — y coordinate 
system. 

The system under consideration is an example of a ma- 
terial whose conductivity is a function of position, i. e. a 
= cr(r), where in the present case r is a two-component 
position vector. In such a material, the macroscopic 
properties can be described by an effective conductiv- 
ity (Je, where (jg is position-independent. If the material 
of interest is additionally anisotropic, as in the present 
case, then cr(r) is a d-dimensional tensor (i. e., a d x d 
matrix, where d is the spatial dimension). In this case, 
(Je remains position-independent but will also become a 
d-dimensional tensor, defined by the relation (see, e. g.. 



(J)=CTe(E). 



(7) 



Here the triangular brackets (....) denote a space- average, 
and J and E are vectors denoting the position-dependent 
current density and electric field, respectively. Thus, for 
example. 



(J> = 1 d''r3{v) 



(8) 



in d dimensions, where V is the (d-dimensional) volume 
of the system. A similar equation holds for (E). 

In general, if cr(r) is a random function of position, <Je 
can only be calculated approximately. However, if cr(r) is 
a non-random function of position, one can often calcu- 
late the needed space-averaged fields exactly, by solving 
the relevant electrostatic equations V-J = 0, VxE = 0, 
J(r) = (j(r)E(r). This is the situation for the case of 
periodic stripes considered in the present paper. 

To get the effective conductivity tensor, we then pro- 
ceed with the following steps. First, we transform the 
conductivity tensors a a and gb into the rotated coordi- 
nate system x' — y' . Next, we calculate the effective con- 
ductivity tensor in the rotated coordinate system. This 
is straightforward by imposing the boundary conditions 
that the component of electric field parallel to the stripes, 
and the component of the current density perpendicular 
to the stripes, should be continuous. Finally, we trans- 
form the resulting conductivity tensor back to the x — y 
coordinate system, in which experiments measure the 
conductivity (or resistivity). 

First, the conductivities and cr'^ in the rotated co- 
ordinates have the following relation to the conductivities 
in the original coordinates a a and <jb, 



(^A 



R-^cfaR 
R ^cfbR-) 



(9) 



where is a 2 x 2 rotation matrix with elements 



R=\ ^"r^. . (10) 

Carrying out this matrix multiplication, we obtain 



cos^ + sin^ ( 



ar sin' + a?^ cos' ( 



= ^ =(af^-af)cos^sin^, (11) 

where i = A or B. 

Now, let us study the components of the effective con- 
ductivity matrix in the rotated system. Note that the 
electric field has components in both the x^- and y^- direc- 
tion, while the corresponding components of the current 
density are 



{J) 



y' 



y X 



T-y y 



{E) 



y' ' 



(12) 
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FIG. 3: [color online] Effective longitudinal resistivities as a 
function of chemical potential /i when the nematic modulation 
is (a) along a diagonal direction and (b) parallel to the ^-axis. 
The chemical potential serves as tuning parameter, while po 
represents the residual resistivity in the isotropic phase. The 
insets display the typical real space patterns of the modulated 
nematic phases. For diagonal modulation (a) both (red) 
and pf^ (black) are larger than in the isotropic phase. On 
the other hand, for parallel modulation (b), is close to 
po, while is larger than po- See the main text for further 
discussion. 



Imposing the boundary conditions mentioned above, we 
obtain the effective conductivity matrix a'^ as follows, 



X y f X y X y \ 



^y'y' 



y'y' 



pay +(l-p)(7^ 
1 

p/ay;y' + {i-p)/ay^y 



(13) 



The final step is to transform the effective conductivity 
tensor back into the original coordinate system, which 
coincides with the crystal lattice. This is done by writing 

(Je= R(j'eR~^, (14) 
where <jg is the matrix whose elements are given in Eq. 



(13). The final result gives the effective conductivity ten- 
sor in the original x — y coordinate system, which will 
depend on the elements of the conductivity tensors of 
the two possible nematic orientations, as well as the pa- 
rameters p and 9. Note that cTg, unlike the conductivity 
tensors of the uniform nematic phase, has nonzero off- 
diagonal elements in the x — y coordinate system. This 
apparently arises because of the asymmetrical geometry 
produced by the stripes, which are at a finite angle to 
the X-axis. 

As we have shown above, the modulation vector q lies 
either diagonal (6> = 7r/4) or parallel {0 = 0) to the crystal 
axes. In both cases, the modulation is determined by 
a single wave vector, so that p=l/2 due to the equal 
population of A and B type. This leads io cf\^ = ^ 



and a' 



X y 



CFj^^ , which simplify the 



above equations to 



1 



( ^X X 



/ X y 



^ A 



^yy 



(15) 



^y y 



2ar 



x'x' 
A 



- (J' 



y y 

A 



x'x' 
A 



a 



y'y' 



For 6 = 7r/4, the effective resistivity tensor using 
^ is then given by 



1 



^yy 



n^y — 



XX 

A 



^A 



^yy\2 



^yy 



2 



XX 

A 



a 



yy\2 

A ) 



XX 

A 



^yy 



Aa^^af (a; 



XX 

A 



^yy 



(16) 



A schematic plot of p^ is displayed in Fig. [s] (a) , where 
the effective resistivities in both directions p^^ and /)|^ 
are higher than the isotropic resistivity pQ. The absolute 
magnitude of p^^ depends on and . For example, 
if cr^ = 0.5(70 and = l.Sao where ctq = Pq^ is the 
conductivity in the isotropic phase, one obtains p^^ = 
pyy = 2.3po- There is a finite p^^ as well due to the 
diagonal stripe orientation as discussed earlier. 

On the other hand, for ^ = as shown in Fig. [s] (b), 
the effective resistivity tensor is given by 
9 XX I ^yy 

pyy ^ ^ o'^y = o fir) 

yy^ He r, xx ^yy ' ^- \^ ' J 

^A ^^A ^A 

Note that the effective resistivity parallel to q, i.e., p^^ 
in this case, is given by the average of two resistivities 
{pyy = py^ + py^ = l/o-^ + l/a^), since the resistiv- 
ities add in series. On the other hand, p^^ = l/cr^^ = 
2/{(j^ + (J%^) is obtained from the average of the con- 
ductivities. It is striking that for the parallel modulated 
nematic state only becomes larger than po, whereas 
pxx Yemmns largely unchanged from the isotropic value. 
This is in contrast to the pure nematic phase, where the 
conductivity along one lattice direction is smaller than 
along the other lattice direction, while the average of 
both is close to the isotropic conductivity, [HI [121 El] 
resulting in one principal resistivity component smaller 
and the other larger than the isotropic resistivity po- 

Discussion and summary - The residual resistivities in 
the spatially modulated nematic state are qualitatively 
distinct from the uniform nematic phase. In the uni- 
form nematic phase, the the Fermi velocities for mo- 
menta in X- and ^/-direction differ (in magnitude) due to 
the deformation of the FS, which then leads to different 
longitudinal resistivities p^^ 7^ pyy. However, in gen- 
eral the resistivity in one direction is higher and that 
in the other direction is lower than the resistivity in the 
isotropic phase po- 

In the modulated nematic phase the resistivities de- 
pend on the angle between the single ordering wave vec- 
tor q and the crystal x-axis, and the strength of the ne- 
maticity given by a^—cr^ in a nematic domain. When q 
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is along the diagonal direction of the square lattice, the 
effective resistivities along x- and y-axes are equal and 
higher than in the isotropic phase. This behavior resem- 
bles the results obtained in Ref. [33j, where the coupling 
between nematic ordering and a specific phonon mode re- 
sults in diagonal nematic domains, which then enhances 
the resistivities p^^ and p^^ . On the other hand, when 
the ordering wave vector is parallel to the crystal axes, 
the direction parallel to q is more resistive, while the di- 
rection perpendicular to q exhibits a resistivity close to 
that of the isotropic phase. 

Our results show a surprising similarity to the longi- 
tudinal, magnetic field tuned resistivities in the bilayer 
ruthenates Sr3Ru207. !2Tl|23j When the magnetic field 
is applied along the c-axis, the longitudinal resistivi- 
ties in a- and b-direction are higher than in the isotropic 
phase, similar to Fig. 3 (a), within a magnetic field win- 
dow, which is bounded by meta-magnetic transitions. As 
the magnetic field is tilted away from the c-axis, the re- 
sistivity parallel to the in-plane field component remains 
more resistive, but the resistivity along the perpendicu- 
lar direction becomes similar to the isotropic resistivity, 
similar to Fig. 3 (b). Thus it is tempting to suggest that 
the magnetic field acts by effectively tuning the angle 0. 

A further support for a modulated nematic phase with 
a diagonal ordering pattern may be gained from consider- 
ing the bare tight binding FS of Sr3Ru207 in the presence 
of a magnetic field as displayed in Ref. [5]. Focussing 
on the 72 FS sheet (which harbors a vHS near the Fermi 
level), one spin component of the 72 sheets opens up near 
momenta (0, dz7r) and (itTr, 0) for a large field along the c- 
axis, giving rise to a cross-shaped hole pocket as shown in 
Fig. 2 (a) in Ref. [5]. The corresponding 2/cF-lines have a 
similar shape as the contour of this 72 pocket and resem- 
ble the 2/ci?-lines displayed in Fig. [2](c) for p = 1.52 but 
rotated by 7r/4. The leading instability for a modulated 
nematic state therefore could involve momenta q lying 
in the diagonal direction, yielding resistivities consistent 
with the transport measurements. When the magnetic 
field is tilted away from the c-axis, aquiring a finite a- or 
b-axis component, the diagonal hot-spots disappear, and 
parallel nematic stripes emerge instead. 

However, the limitation of the present theory also de- 
serves some discussion. As discussed throughout, the 
FS shape is important in determining the direction of 
q. The FS of the bilayer ruthenate is complex (further 
complicated by the presence of an external field) and dis- 
plays more than one vHS near the Fermi level. [M] Thus 
the competition between different q instabilities needs to 
be investigated in addition to the origin of microscopic 
interaction causing such instabilities. Furthermore, the 
present effective medium theory does not contain mag- 
netic field effects such as cyclotron motion. While the 
microscopic origin of a modulated nematic instability is 
still under investigation, the current study provides a way 
to search for such a phase via the transport properties of 



correlated materials. 

In summary, based on effective medium theory we 
studied the electrical resistivity when nematic order is 
spatially modulated in a stripe-like pattern. The modu- 
lation wave vector q is determined by the FS shape near 
the van Hove points. Assuming that the nematic interac- 
tion has only weak momentum dependence near q = 0, 
it was found that the modulation wave vector is either 
parallel or diagonal to the crystal axes. When q is par- 
allel to one of the crystal axes, the resistivity along the 
parallel direction is higher than in the isotropic phase, 
while the resistivity in the perpendicular direction re- 
mains roughly unchanged. On the other hand, when q is 
diagonal to the crystal axes, both p^^ and are larger 
than in the isotropic phase. The present study suggests 
that it is worthwhile to look for modulated nematic in- 
stabilities starting from microscopic Hamiltonians, and 
motivates further experimental exploration in correlated 
materials. 
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